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Abstract 


In quantum detection theory the optimum detection operators 
must commute; admitting simultaneous approximate measurement of 


noncommuting observables cannot yield a lower Bayes cost. The lower 


bounds on mean square errors of parameter estimates predicted by 


the quantum-mechanical Cramer-Rao inequality can also not be 
reduced by such means. 
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Quantum detection and estimation theory has been developed within the 

conventional framework of quantum mechanics, one of the principal tenets of 

which is that only observables associated with commuting operators can be 

[1-3] 

simultaneously measured on the same system. It has been suggested that 

this formulation is too restrictive, that noncommuting operators can be at 

least approximately measured on the same system, and that to include this 

possibility may permit more effective detection, as measured by a lower average 

[4 5] 

Bayes cost. * J We wish to show that no such improvement can be expected. 

The simultaneous measurement of noncommuting observables has been 
treated by Gordon and Louisell . in order to approximately measure certain 
such observables on a quantum-mechanical system S, it is made to interact for 
a time with a second system A, termed the apparatus. It was shown that a 
suitably defined ideal measurement yielding approximate values of the noncom- 
muting observables can be based on the outcome of measurements of commuting 
observables on the apparatus A, or more generally on both S and A. What we 
must therefore do is apply quantum detection theory — with its restriction to 
commuting observables — to the combined system S + A. 

Suppose we are to decide among M hypotheses Hj, H M . Under 


hypothesis Hj the density operator for the combined system at time t is 
S+A 

Pj (t) in the Schrodinger picture. If at an earlier time tQ the density oper- 

S+A , . , . , , ■[ 7 ] 

ator is Pj (tg), the two operators are related by 

p, (t) = U(t, t 0 ) Pj (t 0 ) U (t, t 0 ), (1) 


U(t, t 0 ) = exp 


ft 

-i Hdt'/fi , 

J to 


where H is the Hamiltonian operator for the combined system S+A and fi is 
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Planck's constant h/27r. The operator U is unitary; that is, with U its 
Hermitian adjoint, UU + equals the identity operator 1. 

Let { II j } be a set of commuting projection operators forming an M-fold 
resolution of the identity, 


M 



(3) 


j=l 

On the combination S + A we are to measure these M projection operators at time 
t, and if the k-th yields the value 1, hypothesis is selected as true.^ ^ 

The average cost is then 


C 


M M 



i=l j=l 


Tr[ Pj S+A (t) JI.], 


(4) 


where is the prior probability of hypothesis Hj and is the cost of 
choosing when Hj is true. Let {llj(t)} be the projection operators that 
minimize C when the system S + A is observed at time t; we call these optimum. 
Then by (1) the operators 


ILj(t 0 ) « U + (t, t 0 ) Hj (t) U(t, t 0 ) (5) 

I 

will minimize C when S + A is observed at time tg. Because of the unitarity 
of U(t, tg), the set { n j ( 1 0 ) } also forms an M-fold resolution of the identity 
into commuting projection operators, and the Hj(tg) are optimum at time tg. 
Since the minimization is carried out over all possible M-fold resolutions of 
identity, the minimum Bayes cost must be independent of the observation 

time t. 


Now let us roll time back to an epoch tg before the system S has come 
into contact with the apparatus A. In the Schrbdinger picture this amounts to 
applying the inverse unitary transformation U + (t, tg) to the state vectors of 
the combined system S + A. Because S and A are independent at this time tg, 
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the density operators must now have the factored form 

p^ +A (t 0 ) = p. S (t 0 ) p A (t 0 ), j = 1, it 

Furthermore, as the apparatus A before the interaction has no information 

A 

about which hypothesis is true, p (tg) in (6) must be independent of j. 


The Bayes cost is now 


M M 


C = Tr P (t 0 ) } j } J C i C.j p. (t 0 ) n. . (7) 

i=l j=l -* 

Since S and A are completely uncoupled, and the state of A is independent 


of which hypothesis iL is true, there is nothing to be gained by observing A. 

S A A 

The optimum projection operators n^(tQ) factor as (to) 1 , where 1 is the 

S 

identity operator for the apparatus A, and the set { IT j (tg)) forms an M-fold 

S 

resolution of the identity 1 for the system S, minimizing the Bayes cost 


M M 


Cj C i:j Tr[p b (t 0 ) n ± ]. 


i=l j=l 


A A 


Since Tr[p 1 ] = 1, the minimum value of C is also the minimum value of C 


in (7) and equals the time-independent minimum Bayes cost C n 


The decision 


among the M hypotheses made at time tg is based entirely on the measurement 


of commuting observables on system S. 


Similar considerations apply to estimating the m parameters 0 - 

Q 

( 01 , 02 » • • • » 6 m ) of the density operator p (0) of a quantum-mechanical system 

S. A version of the Cramer-Rao inequality sets lower bounds to mean square 

[3] 

errors of unbiased estimates of 0i, 0 2 ,..., 0 m . Let be an operator 

As 

whose measurement on S yields an unbiased estimate 0 j of the j-th parameter; 

0j must be an eigenvalue of X j . Although in order to be measured simultaneously 
on the same system the operators Xj must commute, the analysis leading to the 
lower bounds given in [3] does not require commutativity of the operators Xj . 


A 



For the class of commuting operators yielding unbiased estimates of the para- 
meters 8 there will exist lower bounds on the mean square errors, and those will 
be greater than or at least equal to the bounds derived in [3] . 

Again including the possibility of measuring noncommuting observables 
cannot lead to lower bounds smaller than those in [3]. In order to measure 
such operators even approximately, a measuring apparatus A must be allowed to 
interact with the system S, and according to Gordon and Louisell's treatment 
of the process, commuting operators will at the end be measured on the combined 
system S + A. In the Schrbdinger picture the density operator p (0, t) 
for S + A will have a time dependence similar to that in (1) . 

Referring to (7) of [3] we see that the symmetrized logarithmic deriva- 
tives (SLD) Lj(t) appropriate for determining the Cramer-Rao lower bounds when 
the measurements of Xj are made at time t are related to those appropriate for 
measurements made at t 0 by 

L-j(t) = U(t , t 0 ) L j ( t Q ) U + (t, t 0 ). (9) 

Then (13) of [3] shows that the matrix A that sets the lower bounds is indepen- 
dent of the time t of observation, again because of the unitarity of the 
operator U(t, to) . 

Once more we move back to an epoch to before the system and the apparatus 

have interacted. The density operator p^ + ^(0, t 0 ) factors as p^(0, to) p^(to), 

A 

where the density operator p (to) of the apparatus A is independent of the 
estimanda 0. The SLD operators for calculating the lower bounds are now the 
solutions of the operator equations 

3 P S (0, t o ) / a e n = | [p s (8, t 0 ) L n s + L n s P s (0, t 0 )], (io) 

A 

and they act only on system S, commuting with p and all other operators on the 
apparatus A. When taking the trace over the states of A to form the elements 
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of the matrix A, the density operator p is replaced by 1, and the lower 

s 

bounds depend only on p (0, to). Thus the apparatus A cannot help estimate the 
parameters 0 of S with smaller mean square errors than the lower bounds cal- 
culated by the quantum-mechanical Cramer-Rao inequality as applied to the 
density operator of system S alone. 

In [3, p. 238] lower bounds were calculated for unbiased estimates m x and 
my of the components of the complex amplitude p = + im y of a simple harmonic 

oscillator, which might represent a mode of the field in an ideal receiver in 
the presence of thermal noise. Those bounds are 

Var \ (N + ~) , Var m y > \ (N + |) , 

where N is the mean number of noise photons. The noncommutativity of the SLD's 
L„ and L„ used to derive these bounds does not invalidate them. It can be 

x y 

shown that if the mode is coupled with an ideal amplifier whose gain is high 
enough to raise the oscillator variables to the classical domain where they 
commute, error variances 

Var m x = Var m y = — (N + 1) 

can be attained [8], It is unknown whether commuting operators can be found 
whose measurement will yield unbiased estimates m x and m y with variances lying 
between — (N + -^) and (N + 1) . 

The measurements we need to make on a quantum-mechanical system S for 
testing hypotheses or estimating parameters will always have to be effected by 
means of an auxiliary apparatus A, and this apparatus, subject to thermal and 
quantum fluctuations of its own, will ordinarily introduce additional random 
uncertainties. Each measurement procedure will have to be analyzed to determine 
what error costs it entails. Detection theory and estimation theory seek 
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lower bounds on these costs, and in doing so they minimize with respect to the 
entire class of possible detection or estimation operators that can be applied 
to the system. The resulting bounds are independent of the time of observation, 
and they cannot be reduced by using any auxiliary apparatus that initially 
posseses no information about the state of the system. 


7 
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